The flexural modes of graphene on a substrate 
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Out of plane vibrations are suppressed in graphene layers placed on a substrate. These modes, 
in suspended samples, are relevant for the understanding of properties such as the resistivity, the 
thermal expansion coefficient, and other. We present a general framework for the study of the 
properties of out of plane modes in graphene on different substrates. We use the model to estimate 
the substrate induced changes in the thermal expansion coefficient, or in the temperature dependence 
of the resistivity. 



Introduction. Since its isolation in 2004 [1 graphene, a 
monolayer of carbon atoms arranged in a honeycomb lat- 
tice, has received great attention due to both its unique 
electrical and mechanical properties [2j. In graphene, the 
carbon atoms display an sp 2 hybridization, with the out 
of plane p z orbitals forming a ir band which is responsi- 
ble for its electrical properties, while the the sp 2 orbitals 
form strong a in plane bonds that govern its mechanical 
properties. It has been verified both experimentally and 
theoretically [3-5] that graphene is the known material 
with the largest in plane elastic constants. 

In freely suspended graphene samples the vibrations 
of the lattice can be classified into in plane and out of 
plane (flexural) modes, with the flexural modes lying at 
lower energies and showing a quadratic dispersion. An- 
harmonic effects at long wavelengths strongly couple in 
plane and flexural modes. The out of plane modes are 
responsible for the significant temperature dependence 
of the electronic resistivity at low temperatures [SHIP]. 
Anharmonic effects in suspended graphene explain its 
negative thermal expansion coefficient, and they play a 
significant role in the thermal conductivity [TT] . 

The interaction between a graphene layer and a sub- 
strate underneath change significantly the properties of 
the out of plane vibrations of the entire system [TTJ [12] . 
Coupling to a substrate also leads to heat transfer be- 
tween the two systems [T3HT5] . 

We present here a simple but generic model which al- 
lows us to study the main properties of a graphene and 
substrate system in the long wavelength limit. We use 
the theory of elasticity, and describe the coupling be- 
tween the graphene layer and the substrate by a simple 
function with a quadratic dependence on the separation. 
As an example we will present results for two of most 
used substrates: silicon dioxide, Si02, and hexagonal 
boron nitride, hBN. We studied the thermal expansion 
of graphene on a substrate and also the effect of the flex- 
ural modes, modified by the coupling to the substrate, 
on the electrical resistivity. 

The model In our model the graphene membrane is 
supported by semi-infinite substrate that occupies the 
half-space z < 0. We will use continuum elasticity theory 
to model the system. We assume that only the out of 



FIG. 1. (a) Pictorial representation of the model used for the 
membrane-substrate coupling, (b) Density plot of the spec- 
tral function A(q, uS) for graphene on Si02 measured in units 
of 7rcJo/(27o)- The straight dashed lines represent from top 
to bottom: uj — c^Q and uj — cr^. The vertical line is at 
q — en/a. The dashed parabola shows the flexural disper- 
sion relation for a free membrane cjf = aq 2 . The solid red 
line, shows the dispersion relation of the hybridized flexural- 
Rayleigh mode, cjfR. 



plane displacements of the membrane couple to the out of 
plane displacements of the substrate as done in Ref. [13] . 
We will assume that the coupling is quadratic and neglect 
anharmonic effects. The action describing the system is 
given by S = ft + ft ubs + ft_ su bs, with 
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d 2 x[^p 2D h 2 -^(\7 2 h) 2 



ftubs = dt d 2 x / dz I -p3DU 2 - -XijkidiUjd k ui ) , 
ft-subs = — / dt I d 2 x (h - u z f . (1) 



dt I d x(h — u z y 

Jz=0 



ft is the quadratic action for a free membrane, where 
h is the out of plane displacement field of the mem- 
brane, p2B is the mass density per unit area and k, is the 
bending rigidity. Form this action one obtains the well 
now quadratic dispersion relation for flexural phonons, 
Wf(q) = aq 2 , with a = \/ft/p2D- For graphene we have 
P2D = 7.6xl0 _7 kg/m 2 and k « l.leV 0. ftubs is the ac- 
tion for a elastic medium occupying the half-space z < 0, 
with mass density per volume p3D, constants given by 



Xijki and displacement field Ui. Finally, 5f_ sll bs describes 
the harmonic coupling between the out of plane displace- 
ment of the membrane and the out of plane displacement 
of the substrate surface, where g is a spring constant per 
unit area. We neglect here the shear between graphene 
and the substrate. Note that, if the substrate is not flat, 
such as Si02 , graphene will have regions that are strongly 
coupled to the substrate and will conform to its config- 
uration, while other regions will become detached from 
it [16] . The detailed analysis of this situation is beyond 
the scope of the present paper, although the analysis pre- 
sented below provides a reasonable approximation to the 
regions which are well attached to the substrate. 

Assuming the dominant interaction is of Van der Waals 
type, g was estimated in Ref. [13] to have a value of 1.82 x 
10 20 J/m 4 for an Si02 substrate. In the absence of more 
rigorous calculations we will use this value both for Si02 
and hBN substrates. The main object of interested from 
which all physically relevant quantities can be obtained 
is the height-height retarded correlation function: 



D(q,u>) 



dtd 2 xQ(t)e 



i(u>t—q-x) 



([h(x,t),h(0,0)]) 



(2) 

where () means thermal and quantum average. Since 
the theory is completely quadratic, it is possible to solve 
it exactly. As a matter of fact, the problem reduces to 
that of two coupled harmonic oscillators. Therefore one 
obtains 



D(q,u) = 



1 



p 2 DU 2 - Kq* - U(q,w) 



(3) 



where II (g, ui) is the correction to the free propagator due 
to coupling to the substrate, 



U(q,uj) 



1 -gA (q,u)' 



(4) 



where Ao(q,uj) is the surface-to-surface substrate propa- 
gator when g = 0: 



A (q,Lj) = - l - \ dtd 2 xQ(t)e* 



(ujt—q-x) 



with both fields evaluated at the surface of the sub- 
strate, z — 0, and where the index o means that g = 0. 
In a quadratic theory, the quantum mechanical prop- 
agator is nothing more than the Green's function for 
the classical equations of motion. Therefore to obtain 
Ao(g, ui) it suffices to study the classical elastic response 
of the substrate to an external pressure at the bound- 
ary z — 0. This analysis was done in Ref. [17] for 
the case of an isotropic elastic medium, such as amor- 
phous Si02, and the result for Ao(#,cj) is given in the 
appendix for completeness. The generalization for the 
case of an uniaxial elastic substrate, such as hBN, will 
be provided somewhere else [18]. The important thing 
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FIG. 2. (a) Plot of the static membrane propagator D(q,0) 
for a Si02 substrate. At high momentum we recover the 1/q 4 
behaviour of the free membrane, but at low momentum the 
behaviour is modified to 1/q due to coupling to the substrate, 
(b) Plot of (/i(0, 0)/i(0, 0)) as a function of temperature for 
Si02 and hBN substrates. As we can see the substrate has 
the effect of stabilizing the membrane. 



to keep in mind, is that a semi-infinite elastic medium 
only supports modes that propagate into the bulk for 
uj > ct#, where c T / L is the transverse/longitudinal veloc- 
ity (replace ct — >• \/C4aJp?>d for the uniaxial case). For 
uj < ct#, the substrate does not support bulk modes, sup- 
porting only a 2D surface Rayleigh mode with dispersion 
given by uj = cr#, where cr is the velocity of the Rayleigh 
surface mode of the substrate. Substrate bulk modes will 
act as a dissipative bath, into which the flexural modes 
can decay, becoming broadened; while the Rayleigh mode 
will hybridize with the flexural mode without any broad- 
ening. Therefore the effect of coupling to the substrate 
on the flexural modes will be twofold: part of the spectral 
weight will hybridize with the Rayleigh mode (which has 
already pointed in Ref. [11]) having a dispersion relation 
approximately given by cjfR ~ c^q and which we will refer 
to as flexural- Rayleigh mode, fR; and part of the spec- 
tral weight will couple to the bulk modes of the substrate 
becoming gapped and broadened, with we will refer to as 
flexural-gapped, fG. This situation is illustrated in Fig. [I] 
where we show a density plot of the spectral function, 
A(q,ui) = —2por>ujImD(q,uj)/7rfl9\, along with the dis- 
persion of the fR mode. The gap and broadening of the 
fG mode are controlled, respectively, by the quantities 
^o = VVP2D and 70 = g/(c L p 3 r>). To obtain the gen- 
eralization to the uniaxial case replace cl — >• y/css/ps^. 
For the used value of g one obtains ujq « 10 meV and 
70 ~ 13 THz for graphene on Si02 and 70 ~ 23 THz 
for graphene on hBN. This values are an overestimation 
comparing with the inverse relaxation times obtained 
from molecular dynamics simulations for acoustic flexu- 
ral phonons on a Si02 substrate, 1/r ~ 0.1 — 1 THz [T2] . 
The hybridization between flexural and Rayleigh modes 
is strongest for q < cr/c^, being more relevant for values 
of the spring constant g such that co>o ~ c r/(2c0- 

Results. Having obtained the propagator, D(q,u), one 
can compute all physically interesting quantities. To 
study the ordering of the membrane on can look for in- 
stance at the same-time height-height correlation func- 



tion: 



(h(x,0)h(0,0)) = -h 



I 



cPqduj 
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2k B T 



lmD(q 



which reduces 
(h(x,0)h(0,Q)) 



in the high temperature limit 

d 2 q 



(6) 
to 



k B Tf-^e iq - x D(q,0), in accor- 



(27T) 



other quantity that is also interesting to look at and that 
is of greater experimental importance is the areal ther- 
cj)pal expansion coefficient of graphene, a a, on a substrate. 
Generalizing the quasi-harmonic treatment of coupling 
between flexural and stretching modes of a membrane 
of Ref. [20] for the case of flexural phonons modified by 
interactions one obtains 



dance with the equipartition theorem. For the quadratic 
dispersion relation of flexural phonons it is known that 
the same point correlation function, (ft(0, 0)h(0, 0)), di- 
verges at any finite temperature, indicating that there 
should be no crystalline order. As expected, the cou- 
pling to a substrate stabilizes the membrane. In Fig. [2] 
we can see a plot of (ft(0, 0)ft(0, 0)) and of D(q,0). It is 
interesting to notice that contrary to what one would 
naively think D(q,0) does not tend to a constant for 
small q but instead goes as 1/q. This is because the 
surface substrate propagator at zero frequency and small 
momentum behaves as A (g, 0) = —l/(Kiq) (with K\ — 
2/>3dCt ( c l ~ c t) / c l f° r an isotropic medium) dominat- 
ing D(q,0). This behaviour will make all the results 
for physical quantities finite, without the need to im- 
pose any infrared cutoff due to anharmonic effects or fi- 
nite sample size. For large distances the height-height 
correlation function now behaves as (h(x,0)h(0,0)) ~ 
ksT '/ , (2ttKix), a result that is independent of g and co- 
incides with the one for (u z (x,Q)u z (Q,0)} at z — 0. An- 
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lm.D(q : u) 



k B T*J (2n) 3H 4sinh 2 
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The obtained value is negative since sgn (ImD(q,uj)) = 
— sgn(cj) and tends to a constant in the high temperature 
limit. Close to room temperature one obtains a value in 
the order of —6 to —7 x 10 -6 K~ a value that is smaller 
in absolute value than the one obtained for a suspended 
membrane [20] as expected. Finally, knowing the form 
of the propagator D(q,u) one can also study the contri- 
bution to the electrical resistivity from flexural phonons 
on supported graphene. Generalizing the results from 
Ref. [8, 9 we compute the resistivity from the known for- 
mula p _1 = \N(€f)v 2 f Tf , where Af(ep) — 2kp/ (ttHvf) 
is graphene density of states at the Fermi level, with vp 
and kp the Fermi velocity and momentum respectively, 
and Tp the transport scattering time. In the quasi-elastic 
approximation, one obtains 



d 2 qd 2 p f duodv 



-h 2 \wk Fi q iP \ (1 - cos 0k Ft k F +q+ P ) d{u ,v)5(c k+q+p - e F )lmD (p, v)lmD(q, u) 



(8) 



where \w k ,q lP \ = q 2 p 2 V 2 {\q + p 
electron-flexural phonon coupling, 

h 2 vl(3~ 



is the squared 
with V 2 (Q) = 



the generalized deforma- 



tion potential [8] that takes into account coupling via 
screened deformation potential and bond stretching vec- 
tor potential [2) [2TU23] , where D(Q) is the screened de- 
formation potential and /3 = — ^ gt a describes the change 
in electron hoping with bond stretching. We have defined 
the function due to thermal occupation factors: 

c ' w ' v ' 4cosh {v) sinh {&) sinh (A) ' 

(9) 
We use a value of Do ~ 25 eV for the unscreened defor- 
mation potential and /? « 2.5 for the bond stretching pa- 
rameter. The electron scattering due to flexural phonons 
is a two phonon process. Therefore there are three contri- 
butions: (i) scattering by two fG modes, (ii) scattering by 
two fR modes, and (iii) a mixed process involving scatter- 
ing by one fR mode and one fG mode. In Fig. [4] we plot 



r 



of the resistivity due to flexural phonons on supported 
graphene on Si02 and hBN, also showing the individual 
contribution of scattering by two fR modes. We see that 
the contribution from fR modes is smaller in hBN than in 
Si02- This is explained by the smaller hybridization with 
the Rayleigh mode in hBN due to a smaller Rayleigh ve- 
locity when compared to Si02- Nevertheless we see that 
the at room temperature the obtained resistivities are of 
order ~ 1Q much smaller than the expected contribution 
from flexural modes in suspended samples [8 and the 
contribution from in plane phonons [9] [24] . 

Conclusions. We have employed a simple model to 
do a complete study of the dynamics of flexural modes 
of a membrane coupled to a substrate, with the aim of 
understanding the role that this modified modes might 
have in the physics of graphene. Since a half-space elastic 
medium supports both 3D bulk modes and a 2D surface 
Rayleigh mode, coupling of the membrane modes to the 
substrate will lead to a splitting of the spectral weight of 
the flexural modes in two parts: one that will hybridize 
with the substrate Rayleigh mode acquiring an almost 




FIG. 3. Areal thermal expansion coefficient as a function of 
temperature for a graphene membrane supported by SiCb and 
hBN substrates, solid line, and its high temperature limits, 
dashed line. It is also show the expected value of the areal 
thermal expansion at high temperature for a free membrane 
imposing a infrared cutoff due to anharmonic effects [20] 



linear dispersion relation for momentum q < CR/a, while 
the other part will become gapped and broadened due 
to coupling to the bulk modes. As expected coupling to 
the substrate leads to a stabilization of the membrane 
and all correlation functions becomes finite, while for a 
free membrane they are known to have infrared diver- 
gences in the harmonic theory. It is worthwhile noticing 
that at high temperature the low momentum behavior 
of the height-height correlation function changes from 
1/q 4 to 1/q. This implies that for large distances the 
height-height function will go as 1/x. We also explored 
the behavior of the areal thermal expansion coefficient 
of graphene on substrate. At room temperature we ob- 
tained a value of the order of —6 to —7 x 10 -6 K _1 . 
Finally we studied the contribution of flexural modes on 
the electrical resistivity of graphene supported by a sub- 
strate. We found that coupling to the membrane strongly 
suppresses the contribution of the flexural phonons, even 
if one takes into account the contribution coming from 
the hybridized flexural- Ray leigh mode. Note that the 
model describes a flat graphene layer on a flat surface. 
This is a good approximation for substrates such as hBN. 
In the case of corrugated substrates, such as Si02, our 
analysis is expected to describe the regions where the 
graphene layer and the substrate are flat and the two 
systems are in close contact. 
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FIG. 4. Computed contribution to the resistivity due to flex- 
ural phonons on supported graphene on Si02 and hBN sub- 
strates for a electronic density n — 10 12 cm -2 . The dotted 
lines show the contribution from the scattering involving two 
fR phonons. 
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C13 





C33 


27.0 


C44 


7.7 


c R (m/s) 


3392 


c R (m/s) 


1835 



TABLE I. Elastic properties for amorphous Si02 and hBN. 
Data for hBN taken from Ref. [25]. 



Substrate surface-to-surface propagator 

The response to an external pressure a z ex at the 
boundary, z = 0, of an elastic medium occupying the 
half-space z < is described by 



u z (q,u) = -A (q,u)ai ext) (q,u), 



(A.10) 



where q is the in plane momentum and u z (q,u) is eval- 
uated at z = 0. The response function, Aq(#, cj), which 
coincides with the quantum mechanical propagator ^ 
for an isotropic substrate is given by [17] : 



Ao(q,u) 



c^byq^uj) 




q 2 + sgn(co>)i0+, 



(A.11) 



with 



S(q,u>) 




— ) - Q 2 + sgn(cj)iO+\ ( — 
ctJ V \ c l 



q 2 + sgn(cj)zO+. 



(A.12) 
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